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RESUMEN 

Favor de proporcionar un resumen en espanol. If you cannot provide a Spanish abstract, the 
editors will do this. We present hydrodynamical models of circumstellar medium (CSM) of long gamma-ray 
burst (GRB) progenitor candidates. These are massive stars that have lost a large amount of mass in the form 
of stellar wind during their evolution. 

There are two possible ways to probe the CSM of long GRB progenitors. 
Firstly, the GRB afterglow consists of synchrotron radiation, emitted when the GRB jet sweeps up the sur- 
rounding medium. Therefore, the lightcurve is directly related to the density profile of the CSM. The density 
can either decrease with the radius squared (as is the case for a freely expanding stellar wind) or be constant 
(as we would expect for shocked wind or the interstellar medium). 

Secondly, material between the GRB and the observer will absorb part of the afterglow radiation, causing 
absorption lines in the afterglow spectrum. In some cases, such absorption lines are blue-shifted relative to the 
source indicating that the material is moving away from the progenitor star. This can be explained in terms of 
wind interactions in the CSM. We can use the CSM of these stars to investigate their prior evolutionary stage. 

ABSTRACT 

We present hydrodynamical models of circumstellar medium (CSM) of long gamma-ray burst (GRB) progenitor 
candidates. These are massive stars that have lost a large amount of mass in the form of stellar wind during 
their evolution. 

There are two possible ways to probe the CSM of long GRB progenitors. 
Firstly, the GRB afterglow consists of synchrotron radiation, emitted when the GRB jet sweeps up the sur- 
rounding medium. Therefore, the lightcurve is directly related to the density profile of the CSM. The density 
can either decrease with the radius squared (as is the case for a freely expanding stellar wind) or be constant 
(as we would expect for shocked wind or the interstellar medium). 

Secondly, material between the GRB and the observer will absorb part of the afterglow radiation, causing 
absorption lines in the afterglow spectrum. In some cases, such absorption lines are blue-shifted relative to the 
source indicating that the material is moving away from the progenitor star. This can be explained in terms of 
wind interactions in the CSM. We can use the CSM of these stars to investigate their prior evolutionary stage. 

Key Words: GAMMA RAYS: BURSTS — STARS: WINDS, OUTFLOWS — STARS: CIRCUMSTEL- 
LAR MATTER — STARS: WOLF-RAYET — LINE: IDENTIFICATION — HYDRODY- 
NAMICS — METHODS: NUMERICAL 



1. INTRODUCTION 

The generally accepted model for lon g gamma- 
ray b ursts (GRBs) is the collapsar model (jWooslev 
19931 ). According to this model a gamma-ray burst 
occurs if a rapidly rotating Wolf-Rayet star (the fi- 
nal evolutionary stage 1 of a massive star) collapses 
toward a black hole. During the collapse, an accre- 
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tion disk is formed around the black hole. The GRB 
is produced by a rclativistic jet that shoots out of 
the pole, driven by energy generated by the accre- 
tion onto the black hole. As the jet moves outward 
it first encounters the outer envelope of the progeni- 
tor star. Beyond this it moves into the circumstellar 
medium that has been formed by the stellar wind 
of the progenitor (see Fig. []}. The jet hits a region 
of free streaming stellar wind, then moves through 
the wind termination shock into the shocked wind re- 
gion and eventually - provided it has sufficient energy 
to pene trate that far - i nto the interstellar medium 
(ISM) (IChevalier et al.l l2004t lEldridee et all l2005t 
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Ramirez-Ruiz et al.|[2005Hvan Marie et al.ll2006allb[ ). 

The GRB jet sweeps up the surrounding medium 
as it expands, accelerating particles to relativistic 
speeds. As a result, synchrotron radiation is emit- 
ted, which wc observe as the GRB afterglow. There 
are two methods to use the afterglow radiation to 
investigate the surrounding medium. 

1.1. Circumstellar density profiles 
The afterglow lightcurvc is directly related to 
the density profile of the surrounding medium. Nu- 
merical models show that the density can cither 
decrease with the radius squared or be constant 
(IChevalier fc Li I |2000(; IPanaitescu fc Kumar 1 12001 , 



20021 : iChevalier et all \2004) . The latter is surpris- 
ing for we would expect to find a free streaming 
wind region close to the star, which would show 
as p (x l/R 2 . The most likely explanation for 
te occurance of a constant density medium is that 
the jet has already passed beyond the free-streaming 
wind region and has entered the shocked wind re- 
ion where the density is indeed nearly constant 



gion where the density is indeed nearly constanl 
(IWiiers Il200ll ; IChevalier et al.ll2004l ; Ivan Marie et al 
2006allbh . However, this requires the wind termina- 



tion shock to be extremely close to the star, which 
is difficult to understand in the case of a Wolf-Rayct 
stars because they tend to have powerful winds. This 
i9s discussed in more detail in \S O 

1.2. GRB afterglow absorption spectra 
Circumstellar matter between the GRB and the 
observer will cause absorption lines in the GRB af- 
terglow spectrum. These lines provide information 
about the composition of the circumstellar medium 
(CSM). Some GRBs show the same line several 
times. These lines are blue-shifted relative to the 
progenitor indicating a system of discrete velocities 
in the CSM. This can be explained by hydrodynam- 
ical interactions between the stellar wind and the 
ISM or between differen t phases of the stellar wind 
van Marie et"ai1 l2005al lbl. l2006bfh We will discuss 



this in more detail in \S [3l 

2. FORMING A CONSTANT DENSITY 
MEDIUM CLOSE TO A GRB 
Most Wolf-Rayet stars have powerful stellar 
winds, which means that the wind termination shock 
lies far (~ 10 pc) from the progenitor star. If this 
is the case, the GRB afterglow can never be gener- 
ated in the shocked wind material, since the GRB 
jet will have ceased to be relativistic long before it 
has reached this distance. If the afterglow is to be 
(partially) generated in the shocked wind, the termi 



Shocked wind 




R2 



Fig. 1. Schematic view of a circumstellar bubble (not to 
scale). The free-streaming stellar wind passes through 
the wind termination shock (Rl) to enter the hot bubble 
of shocked wind material. The high thermal pressure in 
the hot bubble sweeps up a shell (R2), which expands 
into the ISM. 



2.1. Outside influence 

There are several possible causes for the wind 
termination shock to be close to the progenitor 
star, as was discussed in de tail in our earlier pa- 
pers ( van Marie et al. 2006al fbl). Most of these sce- 
narios involve an increase of the 'confining pressure', 
the pressure that restricts the expansion of the cir- 
cumstellar bubble. For example, by increasing ei- 
ther the density or the thermal pressure of the ISM, 
the expansion rate of the circumstellar bubble de- 
creases. Similarly, if the star has a supersonic ve- 
locity relative to the surrounding medium, the ram 
pressure of the stellar motion performs a similar 
function. Unfortunately, the effect of the confin- 
ing pressure is limited. One needs an increase of 
three orders of magnitude to bring the wind termi- 
natio n shock one order of m agnitude closer to the 
~ 2006ah . This means either a 
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natio n shock has to be at < 0.1 pc ([Chevalier et al 
l2004h . 



star 

very high density (~ 10 4 cm -3 ) or a very high tem- 
perature (~ 10 s K) in the ISM. Stellar motion is 
more effective, especially since it can be combined 
with high ISM density. Even so, the star has to move 
quite rapidly to create a bow shock (keeping in mind 
that it may be moving through an H II region or a 
circumstellar bubble where the local sound speed is 
high because of the temperature). 
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Fig. 2. Distribution of GRBs over the redshift-isotropic 
energy space published previously by (|van Marie et al] 
l2006al ). GRBs whose afterglows show a constant den- 
sity profile (circle) tend to have high energy and oc- 
cur at larger redshifts than the afterglows which are 
formed in a free-streaming wind (square). The ones 
marked as triangles are ambiguous, which means that 
either different groups disagree on the nature of its af- 
terglow, or no group has decided on its nature. If the 
symbol is large, several groups have reached the same 
conclusion. Ref e rences f or these afterglow models are : 



Chevalier fe Lil J2000h. iPanaitescu fc Kumar I 
Panaitescu fc Kumar I (|2002h and lChevalier et al.l 



2001) 



2004) 



2.2. Internal influence 

Apart from these outside influences, it is also 
possible to move the wind termination shock closer 
to the star by decreasing the ram pressure of 
the stellar wind. Since weak stellar winds usu- 



ally occur at low metallici ties (jNugis fc Lamers 



2000l : lEldridge fc Vink ll2006h . this would mean that 
GRBs showing a constant density in their afterglow 
lightcurve should occur more often at higher rcdshift. 
Finally, the penetration of the GRB jet into the sur- 
rounding medium depends directly on the energy of 
he GRB. Therefore, a powerful GRB can penetrate 
much deeper into the CSM and bring the wind ter- 
mination shock within reach even at larger distances. 

2.3. Observational test 

Density slopes around several GRBs have 
been determined from the afterglow light curves 
of a number of bursts 



Panaitescu fc Kumar 2001 



1 Chevalier fc Li 
20021 : 



20001: 



Chevalier et al 



2004 ). Their distribution over the redshift-isotropic 



energy space (see Fig. E?t seems to support our idea 
van Marie et al l2006aE The GRBs with a constant 



density profile occur mostly at high redshift and have 
a high isotropic energy. (We use the isotropic en- 
ergy of the burst, since it is irrelevant whether the 
burst has a high absolute energy, or a narrow beam. 
Both cause it to penetrate deeply into the surround- 
ing medium.) 

If the afterglow is generated partially in the free 
streaming wind and partially in the shocked wind, 
as is quite feasible, one would expect to be able to 
observe the transition form one medium to another 
because of t he jump in the local de nsity. However, 
recent work ( Nakar fc Granot 20061 ) shows that the 
transition may not be visible. 

2.4. Alternative explanation 

Another way to constrain the wind bubble is to 
suppose a period of extremely high mass loss before 
the onset of the Wolf-Rayet phase. Such an outburst 
can be observed in the case of rj Carinae, which lost 
a considerable amount of mass ( M > 10~ 2 M^/yr ) 



during the nineteenth century ([Davidson I 119871 ) . 
Whether such a shell would actually constrain the 
Wolf-Rayet wind bubbl e is not certain. Hydro- 



dynamical calculation s (jGarcfa-Segura et al. 1996; 



van Marie et al 2005b) show that a shell, driven by a 
Wolf-Rayet wind, can break quite easily through the 
shell of shocked Red Supergiant wind material sit- 
ting at the wind termination shock. The latter shell 
contains 10-20 M but apparently does not impede 
the Wolf-Rayet wind driven expansion. 

3. BLUE-SHIFTED ABSORPTION LINES IN 
GRB AFTERGLOW SPECTRA 

Some GRBs show a series of absorption lines 
that arc blue-shifted relative to the progenitor; no- 
tably GRB 021004, which shows a t least six absorp- 



tion lines both in C IV and Si IV (jFiore et al.ll2005t 
IStarling et al. 2005 ). These absorption lines show 



that there is a complicated system of discreet ve- 
lo city features between th e GRB and the observer. 
In van Marie et all (l2005al lbl). we showed that these 
absorption features can be explained by stellar wind 
interactions that take place in the CSM of the pro- 
genitor star, u sing a 40 M Q star as described by 
Schaller et ail (|l992T l. The fastest absorption fea- 



tures (~ 3000 km/s) correspond to the wind from 
the Wolf-Rayet star, whereas the slower features 
(~ 150 — 650 km/s) are caused by fragments of 
the shell that the Wolf-Rayet wind swept up into 
the surrounding medium. The presence of these in- 
termediate velocity components places a time con- 
straint on the evolution of the progenitor star. They 
are comparatively short lived (< 5 x 10 4 yrs), since 
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Fig. 3. The column density in [g/cm 2 ] as a function of 
radial velocity for a 40 Mp> s tar during the Wolf-Rayet 
stage (|van Marie et alj|2005al lbh. The zero velocity com- 
ponent is always present, as the interstellar medium it- 
self provides gas that is standing still relative to the star. 
The Wolf-Rayet wind component at 2200 km/s is visi- 
ble during the entire Wolf-Rayet period. However, the 
150. ..650 km/s component, caused by the Wolf-Rayet 
wind driven shell is short-lived. The column density in 
this figure is the average result for 200 radial gridlines. 

they will eventually dissipate into the surrounding 
gas (see Fig. This means that the progenitor 
of GRB 021004 must have had a short Wolf-Rayet 
phase prior to its explosion. 

3.1. Photo-ionization 

The comparatively low ionization states of the 
ions that produce the absorption lines is not so 
easy to explain. If the ions are in the direct line 
of the GRB, one would expect the atoms to be 
at high ionization sta t es out to very lar ge radii 
(jProchaska et all l2006t lLazzati et ail |2006j) . How- 
ever, evidence so far suggests that GR B 021004 had a 



very strong wind, with a high density (jLazzati et al 
2006h . This would put the wind termination shock 
far away from the star. Alternatively, the wind 
might be highly clumped with optically thick clumps 
shielding part of the ma terial from the ionizing 
photons ( Chen et al.ll2006h . A second explanation 
lies in the possibility that the GRB jet is not ho- 
mogeneous. This may mean that the 7-radiation 
of the burst passed through a smaller area than 
the afterglow. Therefore, the material that is 
absorbing the afterglow was not photo-ionized by 
the b urst ( Starling et al. 20051 van der Horst et al. 
2006h . Even so, the wind termination shock could 



not have been very close to the progenitor star, since 
the afterglow itself produces enough high energy 
photons to ioniz e the gas to high ioniz ation states 
at shorter range ( Prochaska et al. 20061) . 



3.2. Alternative explanations 

It was suggested by Mirabal et al. ( 2003h that 
radiative acceleration could account for the pres- 
ence of blue-shifted absorption lines in a GRB af- 
terglow spectrum. However, for the radiative force 
on the particles to be that powerful one would ex- 
pect higher ion ization states than th ose observed for 
GRB 021004 (IStarling et all 120051 ). Nor can this 



explanation account for the large number of dis- 
creet velocity features. Alternatively, the absorp- 
tion features could be interstellar rather than cir- 
cumstcllar in origin. E.g. the intermediate ve- 
locity lines could be caused by an expanding su- 
perbubble. The velocities are rather high, for su- 
perbubbles in local galaxies typically have expan 
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sion velocities of less than 150 km/ s (IHeiles 
Tenorio-Tagle k Bodenheimer 111988 iMartin Ill99 



but certainly not impossible if the progenitor is sit- 
ting inside a starburst region. However, the chances 
of a single GRB progenitor sitting at the center of 
three expanding superbubbles, each with a different 
velocity are not very large. While a supcrbubblc 
shell may account for one of the absorption features 
it is unlikely to account for all of them. Similarly, the 
high velocity components of the absorption spectrum 
can be a galactic wind rather than the stellar wind 
( Prochaska et al. 20061 ). While this would eliminate 
the problem of the low ionization states mentioned 
before, a galactic wind should produce an asymmet- 
ric absorption feature since the acceleration phase of 
the wind would absorb part of the radiation. More- 
over, it is difficult to explain two discreet absorption 
components at high velocities in this fashion. 

3.3. A possible supernova connection 

Recently, three GRBs were detected that showed 
absorption lines at very high blue-shift relative to 
the progenito r : GRBs 050730, 050922 a nd 060418 



(|D'Elia H20061 : IPiranomonte et al. I [2006al |bl) . If the 



lines observed in these afterglows are indeed pro- 
duced in the GRB progenitor host galaxy, the rel- 
ative velocity of the gas would be of the order 
of 10 4 ...10 5 km/s. Such velocities would be very 
difficult to explain in terms of Wolf-Rayet wind 
velocities, which are usually at least an order of 
magnitude lower. This may indicate, that these 
lines are just foreground noise, occurring in a dif- 
ferent galaxy with a lower redshift than the GRB 
host galaxy. However, there is another possible 
explanation. Velocities such as these do occa- 
sionally occur in supernova explosions, as observed 
in blue-shifted absorption li nes in the spectra o f 
Type II supernova SN2005cs (|Pastorello et al.ll2006f) 
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and Type lb supe rnova SN1987M (lElhamdi etal 



20041) 



See also 



Chevalier fc Fransson 



Mazner fc McKee 
(2006f) 
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For the supernova 
ejecta to appear in absorption in the GRB afterglow, 
the supernova has to happen well in advance of the 
GRB. This is not part of the normal collapsar model, 
but th e possibility was suggested bv lVietri fc Stella 



(|1998l ). who proposed a 'supranova' model, where 
a supermassive neutron star falls back into a black 
hole after the original supernova. The time interval 
between the two events could be as large as several 
years. If some of the high velocity absorption lines 
observed in GRBs 050730, 050922 and 060418 are 
indeed caused by the 'supranova ejecta' this would 
be a considerable step forward in our understanding 
of GRB formation. 

4. DISCUSSION 
The circumstellar medium of a massive star car- 
ries the fingerprints of the previous stages of the 
stellar evolution. It provides an excellent oppor- 
tunity to investigate the evolutionary past of the 
central star. We have presented two methods for 
analyzing the CSM of GRB progenitors. Through 
blue-shifted absorption lines in the GRB afterglow 
spectra we can not only study the composition of 
the CSM, but also the hydrodynamical interactions 
that took place, which in turn give us clues about 
the wind parameters of the previous evolutionary 
phases of the progenitor star. The density profiles 
of the CSM provide us with information on the bal- 
ance between the ram pressure of the stellar wind 
and the confining pressure of the ISM. There is no 
reason, why these tools can not be used in other sit- 
uations as well. E.g., on supernovae and Wolf-Rayct 
stars. If the CSM of these objects can be matched 
to the CSM around GRBs, it will be a major step 
forward in our understanding of the evolution of 
GRB progenitor stars. In the future we hope to 
provide similar models of the circumstellar medium 
around rapidly rotating sta rs ( Yoon fc Langer 1120051: 



Yoon. Langer fc Norman 20061 : Wooslev fc Heger 



20061) which arc thought to be GRB progenitors 
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